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ABSTRACT 


Practical! zing  quadratic  optimal  control  algorithm!  were  used  to  design  load 
relief  systems  for  the  C-5A,  a  large  flexible  aircraft.  The  predicted  rms 
stresses  at  the  wing  root  were  reduced  by  more  than  40  percent.  Handling 
qualities  or  stability  were  not  compromised.  The  control  is  realized  with  a 
gyro  and  three  accelerometers  affecting  ailerons  and  elevator  —  two  accel¬ 
erometers  more  than  an  existing  stability  augmentation  system.  The  quadratic 
performance  index  is  defined  to  enforce  good  handling  qualities  and  to  limit 
the  control  system  bandwidth. 
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SECTION  I 

INTRODUCTION  AND  SUMMARY 


Quadratic  methodology  was  applied  to  the  design  of  a  pitch-axis  load- 
relief  control  system  for  the  C-5A  aircraft.  Predicted  rms  stresr.es  due  to 
wind  gusts  at  the  wing  root  were  reduced  by  50  percent,  while  stress  rates 
were  reduced  by  31  percent.  This  was  done  without  serious  degradation  of 
the  handling  quantities  or  the  stability  of  the  aircraft.  Similar  reductions  in 
peak  and  steady-state  stresses  due  to  maneuvers  were  realized.  Symmetric 
ailerons  and  the  inboard  elevator  are  driven  by  control  signals  from  acceler¬ 
ometers  and  a  gyro.  One  accelerometer  and  the  gyro  already  exist  in  the 
stability  augmentation  system  (SAS).  The  additional  load  improvement  con¬ 
trol  system  (LICS)  acceleiometers  would  be  placed  on  outer  wing  panels. 

The  effectiveness  of  active  control  for  effecting  load  improvement  is 
summarized  in  Table  I.  Results  are  given  for  two  systems.  System  I  uses 
ailerons  to  reduce  the  loads,  whereas  System  n  uses  spoilers  in  addition  to 


Table  I.  Wing  Root  Stress  Relief  Summary 


Ratio  of  Controlled  Aircraft  to  Free  Aircraft 

Parameter 

Syatem  I;  Atlorona,  One 
Extra  Sonaor  Sot 

Sjratem  U:  Allarona  +  Spotlora, 
Toro  Extra  Sonaor  Setac 

Wing  root  itron* 

0. 50 

0.23 

Wing  root  atroao  rate* 

0.69 

0.59 

Peak  manouroiiag 
atroao**" 

0.5T 

0.36 

Steady -at  at* 
maneuvering 

atroaa** 0 

0.57 

0.34 

*For  Lockheed  night  condition  ST  <w  ■  593, 154  lb:  M  ■  0.53S;  h  ■  10,000  ft) 
using  •  itx- flexure -mode  rep  re  notation. 

bTor  *1. 5  Incremental  g  cuing  25  dog  of  tap  afloron.  Tbo  aileron*  would  hit 
tho  down  (top*  o f 15  dog  at  -0. 9  incromantal  load  factor. 

eFor  10  dig  of  up  apotlor  at  +1.5  Incromantal  g. 
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the  ailerons.  Table  I  shows  that  rms  stress  and  stress  rate  can  be  reduced  by 
50  percent  and  31  percent,  respectively,  using  System  1.  Peak  and  steady- 
state  maneuvering  stresses  can  be  reduced  by  43  percent  of  their  nominally 
attained  values  using  System  I.  System  I  requires  one  additional  set  of  sen¬ 
sors  to  achieve  the  gust  relief  improvements  cited.  In  the  study  conducted, 
virtually  no  improvement  was  obtained  using  additional  sensors. 

By  using  spoilers,  performance  could  be  Improved  as  indicated  for  System 
n  in  Table  I.  For  System  n,  use  of  a  pilot-operated  switch  is  recommended  to 
activate  the  spoiler  portion  of  the  system.  In  normal  use  the  spoilers  would 
not  be  deflected;  the  performance  would  be  that  of  System  I.  In  rough  air, 
the  spoilers  would  be  activated  to  achieve  the  results  noted  for  System  n.  The 
spoilers  would  be  biased  at  10  degrees  in  this  mode  of  operation. 


•*» 
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SECTION  n 

P  ROBLEM  FORMULATION 


DESIGN  PROCEDURE 

The  design  approach  used  is  based  on  quadratic  optimal  control  theory  -- 
optimal  with  respect  to  a  quadratic  performance  index  subject  to  practical  con* 
straints.  The  background  for  this  methodology  is  given  in  References  1  and  2. 

The  design  is  accomplished  by  minimizing  a  quadratic  performance  Index 
which  weights  mean  square  stresses,  stress  rates,  model- following  errors, 
control  surface  rates  and  control  surface  deflections.  Simple  compensation 
filters  were  included  in  the  measurement  constraints  to  improved  bending¬ 
mode  damping  and  handling  qualities.  The  model-following  errors  were 
weighted  to  enforce  good  handling  qualities. 

The  quadratic  optimal  design  technique  requires  that  the  aircraft  be 
modeled  as  a  linear  time-invariant  plant  representing  a  single  flight  condition: 


• 

X 

•  Fx  +  GjU  +  Gjti 

(1) 

r 

«  Hx  +  Du 

(2) 

y 

-  Mx 

(3) 

where 

x  ■  State  vector  (including  rigid-body  states,  actuator  and  servo 
states,  flexure-mode  states,  sensor  states,  model-following 
states,  and  wind  states) 


S 


u  ■  Control  input  vector 
r)  ■  Unit -variance  white  nolae  vector 
r  ■  Response  vector 
y  ■  Measurement  vector 

CONTROL  LAW 

The  control  law  is  constrained  to  be  of  the  feedback  gain  form 

u  «  K*y  (4) 

where  the  asterisk  is  used  to  differentiate  this  matrix  K*  from  the  optimal 
full-state  foedback  form 

u  «  Kx  (5) 


PERFORMANCE  INDEX 

The  performance  index  is  defined  to  be 

J  -  E{Tr[Q  r  rT]  }  (6) 

THE  PROBLEM 

Using  Equations  (1)  through  (6),  the  problem  reduces  to  minimising  the 
performance  index  J  with  respect  to  the  gains  matrix  K*  subject  to  the  con¬ 
straint  of  system  stability  and  Equations  (1),  (2),  and  43). 
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There  are  two  parte  to  this  problem.  Thia  firat  part  involves  determining 
the  full-state  feedback  C  Equation  (5)]  .  The  second  part  Involves  constraining 
the  feedback  [Equation  (4)]  and  is  also  known  as  the  fixed-form  optimal  con¬ 
trol  problem.  The  first  part  of  the  problem  is  discussed  in  Section  El.  The 
simplified  control  law  is  discussed  in  Section  IV. 

AIRCRAFT  MODEL 

The  design  approach  was  applied  to  the  pitch  axis  of  the  C-5A  aircraft 
at  one  flight  condition  (FC-37)  which  is  a  low-altitude  cruise  condition  with 
the  following  parameters: 

•  Gross  weight  -  593,154  lb  (50l  fuel,  501(  cargo) 

•  Mach  number  -  0. 533 

•  Altitude  -  10,000  ft 

•  Dynamic  pressure  -  292  psf 

•  True  airspeed  -  577  fps 

•  Center-of-gravity  location  -  31<  MAC 

The  general  arrangement  of  the  C-5A  is  shown  in  Figure  1.  The  spoilers 
being  considered  are  the  outboard  spoilers  (which  consist  of  three  panels). 
Figure  2  shows  the  locations  of  the  sensors  that  will  be  considered  in  the  con¬ 
trol  synthesis. 

The  procedure  described  in  Reference  3  was  used  to  reduce  a  C-5A  model 
to  the  form  of  Equation  (1).  The  procedure  included  the  quasi-elastic  effects 
of  bending  modes  neglected,  as  only  the  three  most  prominent  modes  were 
included  in  the  design.  The  final  designs  were  checked  using  a  six-mode 
model. 


The  Lockheed  data  (listed  in  Ref.  4)  for  FC-37  were  processed  by  the 
computer  programs  developed  for  a  com  '-configured  vehicle  (CCV)  (Ref.  3). 
Satisfactory  agreement  with  load  allevla  a  and  mode  stabilization  (LAMS) 
results  (Ref.  5)  was  not  achieved.  Two  errors  in  the  processing  program 
(Table  II- 13  of  Vol.  II  of  Ref.  3)  were  found.  After  making  these  corrections, 
the  results  shown  in  Table  n  were  obtained.  Agreement  is  now  considered  to 
be  sufficient  for  the  intended  purposes. 

The  notation  used  for  column  headings  in  Table  II  (e.g. ,  LAMS-6)  refers 
to  the  data  source  [LAMS  (Ref.  5)  or  this  contract]  and  the  number  of  bending 
modes  used.  The  number  (shown  in  row  2)  is  used  to  convert  0  from  rad / 
sec  to  in/ sec.  It  is  equal  to  0.6066  x  10"^. 


The  state  vector,  x,  for  optimal  state  feedback  designs  was  (17  states): 
xT  ■  (w*  Of*}*  ’Ig*  ’Ig*  ^i*  ^g*  ^3*  fia*  °e  »  Pi*  Pg*  Ps»  ^ 

•V  <V  p6-  V 

where  w  is  the  rigid-body  vertical  velocity  (in.  /sec),  d/n^  if  a  normalized 
rigid-body  pitch  rate,  is  the  first  bending-mode  coordinate,  is  the  sixth 
bending-mode  coordinate,  rig  is  the  third-mode  coordinate,  0a  is  symmetric 
aileron  displacement,  de  is  the  inboard  elevator  displacement,  Pj  through  Pg 
are  wind  distribution  and  lift  growth  states,  and  w^  is  the  vertical  gust  velocity. 

The  control  input  vector,  u,  contains  the  inputs  to  the  actuator  models 
(two  controls): 
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Table  n.  Comparleon  of  RMS  Responses  Due  to  1-FPS 
RMS  Oust 


P*ramit«r 

LAMS-a 

LAMS- 18 

LICS-S 

LICS-6 

w 

(in/  lie) 

10.43 

... 

11.08 

11.03 

i/nj 

(in/  sic) 

l.4t 

... 

1.13 

1.133 

(in.) 

0.431 

... 

0.840 

0.880 

(in/ lie) 

1.61 

... 

1.033 

1.04  1 

(in.) 

0. 0333 

... 

0.0331 

0.0330 

’’s 

(In/ lie) 

o.sn 

... 

0.314 

0.330 

n6 

(in.) 

0.01363 

... 

0.0086 

0. 00066 

<6 

(in/ lie) 

0.1463 

.... 

0.0411 

0. 0486 

•l 

(pi!) 

136.6 

S3*  4 

176.0 

170.0 

*1 

(pal/ lie) 

833.4 

1040.0 

871.0 

864.0 

■j 

(pit) 

143.8 

343.1 

163.8 

181.0 

•a 

(pil/aic) 

863.7 

1876.0 

841.0 

870.0 

No.  of  modoa 

6 

16 

3 

6 

Modn  mid 

1. 3,4.6. 

7.11 

AH 

1.3.6 

1.3. 3. 4.8.8 

Oust  pentt  ration 

No 

Tm 

Yia 

Yia 

Wagner  dynamic* 

Conatant 

(unity) 

Second- 

ordir 

Conatant 

(unity) 

Constant 

(unity) 

Modi  approxi¬ 
mation 

Truncated 

Com  pi  it  id 

Steady- 

atate 

ritalnad 

Stiad-Btati 

retained 

Scali 

Itnfth  (ft) 

1000 

1000 

1780 

1780 

The  complete  transfer  function  for  the  elevator  and  aileron  is  third’ 
order  said  was  used  for  evaluation  only: 


(9) 
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RESPONSE  SELECTION 


The  objective  of  the  LICS  design  was  to  significantly  reduce  wing  root 
stress  and  stress  rate  using  active  control  without  degrading  the  existing 
handling  qualities.  This  reduction  was  to  consider  both  maneuvering  relief 
as  well  as  gust-induced  stresses.  Therefore,  the  response  vector,  r,  mini¬ 
mized  was 

rT  “  a2'  *1*  *2'  8a*  V  6a'  8e*  ®/n2)  (10) 

where  s^  is  the  stress  at  the  wing  root,  Sg  is  the  stress  at  a  mid-wing  panel, 
and  9/  Og  is  a  model-following  error.  The  model -following  error  represents 
the  difference  between  the  pitching  moment  equation  of  the  aircraft  and  a 
rigid  model  with  desired  handling  qualities.  That  is, 

J/ng  -  <F2  -  P2m)  x  (11) 

where  Fg  Is  the  second  row  a!  the  stability  matrix  F  of  the  aircraft  and  Fgm 
is  the  second  row  of  the  stability  matrix  of  the  model.  This  row  is  given  by 


F2m  *  (F21m'  F22m*  0 . °’  F2.11,*"  F2,17) 


(12) 


Here,  E^lm  and  Fggm  are  selected  so  that  the  desired  short-period 
frequency  and  damping  ratio  are  computed  from 


9 

M  *  *  P  F  •  P  P 

usp  rllr22m  r21mr  12 


(13) 


2  ^sp  wsp  *  ~F11  ’  F22m 


(14) 


where  Fjj  and  F^2  are  elements  of  the  first  row  of  F.  This  is  because  the 
z  force  (row  1  of  the  model)  is  taken  to  be 


lm 


(F11'F12'0,  •  •  *  '°'F1, 11,  *  *  •  •  Fl,  17) 


(15) 
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The  wind  coefficients  of  the  model  are  taken  to  be  the  same  as  the  air¬ 
craft  itself.  The  wind  coefficients  of  the  model  in  Equation  (12)  should 
probably  be  adjusted  to  agree  with  an  adjustment  in  F2jm  the  coefficient  of 
w.  However,  in  this  design,  F21m  and  F22m  are  the  8ame  aB  F21  and  F22* 
Thus,  the  model-following  error  tends  only  to  decouple  the  bending  modes 
from  the  rigid  body. 

The  matrices  of  Equation  (1),  (2),  and  (3)  corresponding  to  this  flight 
condition  are  tabulated  in  the  Appendix. 


SECTION  III 

OPTIMAL  STATE  FEEDBACK  CONTROL 


Optimal  full-state  feedback  designs  were  obtained  for  various  quadratic 
weighting  selections  on  the  responses  of  Equation  (10).  Not  all  responses 
were  weighted}  some  were  only  monitored  to  ensure  that  they  were  not  com¬ 
promised.  This  was  especially  true  with  s2  and  s^,  the  mid- wing  stress  and 
stress  rate. 

The  quadratic  weights  selected  versus  chronological  designs  are  shown 
in  Figure  3.  Zero  weights  are  not  shown.  The  weights  on  i  and  0/n«  are 

•  4  *  * 

not  shown,  since  they  were  held  constant  at  20  and  10  .  respectively.  These 

weights,  as  well  as  the  initial  weights  (iteration  1).  were  determined  in 
previous  designs  with  erroneous  stress  equations.  The  weights  on  6  and 
?/ng  were  iteratively  selected  to  constrain  the  rigid-body  frequency  and 
damping  ratio.  This  was  done  with  relatively  little  effect  on  stress  control. 
Figures  4.  5.  und  6  present  the  root  mean  square  (rms)  responses  and  rigid- 
body  frequencies  and  damping  ratios  corresponding  to  the  weights  of  Figure  3. 
The  dashed  lines  represent  responses  which  are  not  weighted.  Table  III  lists 
the  corresponding  frequencies,  damping  ratios,  and  actuator  root  locations  of 
14  optimal  (full- state  feedback)  control  designs.  The  rationale  used  in  the 
designs  will  now  be  discussed. 

On  iteration  1,  the  aileron  actuator  root  was  large  (caused  by  a  signifi¬ 
cant  negative  feedback  of  to  uq),  Thus,  the  aileron  rate  weight,  ,  was 
Increased  for  iteration  2.  There  was  a  corresponding  increase  in  stress  s^. 

In  an  attempt  to  simultaneously  decrease  0  feedback  and  increase  aileron 
effectiveness,  the  weight  on  aileron  displacement,  Qg  ,  was  removed  for 
iteration  3.  This  had  very  little  effect.  Thus,  for  iteration  4,  this  weight 
was  reinstated  and  the  aileron  rate  weight,  ,  was  increased  further. 

The  aileron  feedback  became  positive  and  the  stress  increased;  so,  on  the 
fifth  itr  ration,  the  weight  was  cut  in  half. 
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STRESS  (PSD  AMD  STRESS  RATE  (PSI/SEC) 


Figure  3.  Chronological  Weighting  of  Responses  for  Optimal 
State  Feedback  Control  Design 


NOTE:  ITERATION  ZERO  REFERS  TO  FREE  AIRCRAFT 


Figure  4.  Chronological  Stress  Responses  for  Op  imal  State 
Feedback  Control  Design 
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NOTE*  ITERATION  ZERO  REFERS  TO  FREE  AIRCRAFT 


Figure  5.  Chronological  Model -Following  Responses  for  Optimal 
State  Feedback  Control  Design 


Figure  6.  Aileron  and  Elevator  Responses  for  Optimal 
State  Feedback  Control  Design 


The  iteration  3  controller  was  quite  acceptable,  if  the  5&  feedback  could 
be  successfully  eliminated  later  during  practicalization  designs.  Stress  was 
reduced  67  percent  on  that  iteration,  and  stress  rate  SO.  5  percent.  The 
iteration  5  controller  is  a  compromise,  with  a  49.5-percent  reduction  in 
stress  and  a  42. 3-percent  reduction  in  stress  rate. 

In  an  attempt  to  reduce  stress  rate  even  more,  the  stress  rate  weight, 

,  was  increased  and  Qg  was  removed  for  iteration  6.  Both  stress  and 
stress  rate  decreased,  with  a  significant  increase  in  bending- mode  damping, 
although  the  model- following  error,  short-period  damping  ratio,  and  0 

© 

feedback  to  ug  increased  somewhat.  This  controller  was  also  considered 
a  candidate  for  practical  design;  however,  from  past  experience,  it  was 
expected  to  lead  to  difficulty.  Generally,  practical  measurements  cannot 
produce  this  much  bending-mode  damping  without  excessive  filtering.  This 
is  especially  true  when  slow  actuators  are  used. 

Iterations  7  through  11  were  further  attempts  to  reduce  stress.  Finally, 
by  deemphasizing  reductions  in  stress  rate,  stress  could  be  reduced  much 
more  within  the  constraints  of  the  actuator  bandwidths.  The  stress  rate 
weight,  ,  was  removed  by  iteration  12.  However,  the  0R  feedback  to  ug 
was  still  high.  Iterations  13  and  14  brought  this  feedback  to  within  reason 
for  practical  design.  The  iteration  14  controller  was  also  a  candidate  for 
practical  design.  With  this  controller,  stress  was  reduced  67. 8  percent  and 
stress  rate  36  percent,  with  small  6g  to  ug  feedback. 


SECTION  IV 

SIMPLIFIED  FEEDBACK  CONTROL 
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This  section  summarises  the  LICS  control  design  effort.  The  first  set 
of  controllers  used  ailerons  and  elevator  and  was  based  on  a  three-mode 
model.  Both  full-state  and  measurement  feedback  were  evaluated  (Table  IV). 
Three  of  the  14  full-state  feedback  controllers  were  successfully  simplified. 


The  second  set  of  controllers  also  used  aileron  and  elevator  but  was 
based  on  a  six-mode  model.  .Both  full-state  and  measurement  feedback  were 
evaluated  (Table  V). 


The  last  set  of  controllers  used  spoilers  and  elevator  with  a  three-mode 
model.  Only  full-state  control  was  considered  (Table  VI). 


With  the  iteration  3,  5,  8.  and  14  controllers  as  baselines,  practical 
designs  were  attempted  using  their  respective  quadratic  weights  and  their 
optimal  gains  as  starting  points.  The  practical  designs  were  attempted  with 
different  measurement  complements.  The  candidate  measurements  are  listed 
in  Table  VII. 


For  most  of  the  practical  designs,  the  state  vector, 
states: 


*T.  included  18 


|w,  S/nj#  ^1*  ^8*  ^3*  ^1*  t*3»  ®a*  ®e 

pr  p2*  p3*  p4*  p5»  p6'  wg) 


(16) 
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Table  IV.  Practical  Design  Results  Summary 
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Table  V.  Design  Results  Summary  with  Six-Mode  Model 
towlnd  ■  5,2  fP»> 


Paramatar 

Controller 

Aircraft 

1 

a 

4 

5 

a,  (10s  pat) 

0. 030 

0. 483 

0.311 

0.  470 

0. 581 

0.467 

a,  (I0S  pal) 

0.940 

0.380 

0.943 

0.873 

0. 330 

0. 887 

a,  (10s  pai/aae) 

3,828 

1. 888 

1, 887 

3.381 

3.383 

3.338 

■j  (10S  pai/aac) 

1.93T 

3.487 

3.  579 

3. 137 

3,938 

3. 010 

wtoU,  ((pa) 

43. 78 

43.10 

43.00 

39. 43 

38.78 

39.83 

ilt^  (ln/aac) 

5. 85 

5.80 

5. 11 

4. 55 

5. 15 

4.  35 

ik  (rad/aac) 

0 

0.0488 

0.0784 

0.0871 

0.0438 

0.0088 

flt  (rad/aac) 

0 

0. 0150 

0.0183 

0.0184 

0. 0103 

0. 0154 

8#  (rad) 

0 

0. 0188 

0.0381 

0.  7355 

0. 01031 

0.0388 

8.  (rad) 

0 

0. 0045 

0.0058 

0.0049 

0.0033 

0.0047 

Actuator  modal 

--- 

1* 

1* 

Sb 

3k 

Sb 

Peedbacka 

Nona 

All 

atataa 

All 

atataa 

Same  aa 

sc* 

Tab  la  IV 

Same  aa 

5C  rf 

Table  IV 

Same  aa 

5C  rf 

Table  IV 

*Flret- order  modal  [Rofor  to  Equation  <»]. 
^bird-order  modal  CRafar  to  Equation  (•)]. 


Table  VI.  Spoiler  Effectiveness  for  Gust  Relief 
wind  *  5  2  'l"’* 


Parameter 

Free 

Aircraft 

State 
Control  1 

State 
Control  3 

a,  (10*  pat) 

0.185 

0.037 

0.084 

aj  (10S  pat) 

0.951 

0.321 

0.511 

<1  (10s  pat/aac) 

3.973 

0.388 

0.885 

<2  (10s  pai/aac) 

4.390 

0.944 

1.407 

•total  ^ 

37.84 

43.94 

43.37 

t/nj  (ln/aac) 

3. 81 

7.38 

8.83 

0,  (rad/aac) 

0 

0.0058 

0.054 

0#  (rad/aac) 

0 

0.0198 

0.0148  | 

«,<rad) 

0 

0.0034 

0.0188 

®,  (rad) 

0 

0.0059 

0.0088 

•(  (rad) 

0 

0. 039 

0.0345 

(rad/aac) 

0 

0.094 

0.0598 

"Three-mode  rap raaan tattoo  with  ft  ret -order  actuatora, 
daHactad  apotlar  ajatam,  Mato  control. 
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Table  VII.  Candidate  Measurement  for  Practical  Designs 


Measurement 
Mum  bar 

Description 

1 

Accelerometer  locatad  at  fuselage  panel  4  (forward) 

3 

Difference  between  accelrrometer  located  at  fuselage 
panel  4  and  average  of  accelerometers  located  at  wting 
panels  21  'tip) 

3 

Difference  between  ac  -rierometer  located  at  fuselage 
panel  4  and  average  of  accelerometers  located  at  wtng 
panels  IS  (mid-wing) 

4 

Combination  of  laggsd  and  hlghpassed  rate  gyro  located 
at  fuselage  panel  4  (forward)  with  frequency  cutoff  at 

0.  88  rad/ sec  (-&*) 

5 

Rate  gyro  located  at  fuselage  panel  34  (aft) 

6 

Lagged  measurement  3  with  frequency  cutoff  at  3  rad/ sec 

which  is  the  same  as  Equation  (7)  with  the  addition  of  6 the  lagged  fuselage 
panel  4  rate  gyro  output  used  to  construct  measurement  4  of  Table  VII.  In 
the  designs  where  the  lagged  acceleration  measurement  was  used,  the  state 
vector  included  19  states: 


xT  ■  (w,  9/nj,  i^,  lig.  tir  rig,  rig,  «a,  «e.  9l, 
aL*  Pi*  Pg’  3*  P4*  P5 »  Pg» 
where  aL  is  the  lagged  acceleration  measurement. 


(17) 


The  design  procedure  described  in  References  1  and  2  was  used  to 
realize  the  practical  designs.  Using  this  procedure,  the  measurement  gains 
are  written  as  a  function  of  a  scaler  parameter,  X,  such  that 

K*(X)  «  K1  (X)  +  XK2  0  s  X  *  1  (18) 
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The  Incrementing  parameter,  X,  is  equal  to  1  for  the  optimal  state  feedback 
controller,  and  \  is  equal  to  zero  for  the  optimal  measurement  feedback 
controller. 

The  starting  point  (X  *  1)  is  found  by  using  the  optimal  state  feedback 
gains  and  the  measurement  matrix  (augmented  with  direct  measurements  of 
states  not  necessarily  measurable  so  M"1  exists): 

K*(l)  *  KM"1  (19) 

The  measurement  constraints  are  applied  gradually  by  stepping  X  to  zero. 

The  matrix  K^O)  is  the  fixed-form  solution  and  has  the  gain  structure 
desired. 

This  procedure  of  "backing  off"  from  the  state  feedback  controller  is 
illustrated  in  Figure  7  for  the  design  of  controller  14C  (defined  below).  The 
same  quadratic  performance  index  was  minimized  while  the  measurement 
constraints  were  gradually  applied. 

The  cost  J  is  minimized  at  each  increment  of  X,  with  respect  to  the  gains 
on  the  measurements.  The  gains  are  predicted  for  each  increment.  The 
minimization  of  J  is  the  correction.  The  differences  between  the  gradient 
norms  before  and  after  the  corrections  are  also  shown  in  Figure  7. 

The  practical  designs  are  summarized  in  Table  IV,  except  for  the 
practicalization  of  the  iteration  6  controller.  In  an  attempt  to  practicalize 
this  controller  using  measurements  1,  2,  3,  and  4,  the  rigid- body  ar.d  first 
bending- mode  frequencies  and  damping  ratios  became  unacceptable,  as  their 
roots  appeared  to  merge  with  actuator  roots.  The  successful  practical  designs 
will  now  be  discussed. 
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Figure  7.  Incremental  Gradient  for  C-5A  Design 

Successful  designs  were  achieved  using  Iteration  3,  5,  and  14  controllers 
as  baselines.  Cohtroller  BA  had  one  more  accelerometer  measurement  than 
SB.  This  accelerometer  gained  very  little  stress  and  stress  rate  reduction 
at  the  expense  of  overcontrolling  the  rigid  body.  It  slightly  improved  mode  3 
damping,  with  less  damping  on  modes  1  and  6.  On  controller  5C,  aileron 
feedbacks  were  permitted.  This  defined  a  prefilter  for  the  actuator.  There 
was  an  improvement  over  SB  in  stress  rate  and  stress  Sj  (not  weighted)  at 
the  expense  of  bending-mode  damping.  The  aft  rate  gyro  was  included  in 


controller  5D,  with  some  Improvement  over  SB  in  stress,  stress  rate,  and 
mode  damping. 


Controller?  3A  and  3B  are  baaed  on  the  iteration  3  controller.  On  con* 
troller  3A,  aileron  feedbacks  were  allowed,  with  significant  improvements 
in  stress  and  stress  rate  at  the  expense  of  bending- mode  damping. 

Controller  14  practical  controllers  14A,  14B.  and  14C  produced  lower 
stress  Sj  levels  than  the  others.  Stress  rate  levels  were  higher.  On  con¬ 
troller  14C,  measurement  6  (lagged  acceleration)  increased  bending-mode 
damping  and  lowered  the  short-period  frequency.  The  short-period  damping 
ratio  was  somewhat  higher,  however. 

Controllers  5D,  3B,  and  14C  are  the  most  desirable  controllers  from 
the  bending- mode  damping  viewpoint,  which  could  be  important  for  stability 
margins.  Controller  3B  is  the  least  complex  of  the  three;  however,  14  C  is 
only  different  by  3  lag  network.  Controller  5D  uses  an  extra  sensor. 

Table  V  summarizes  the  performance  of  two  state-feedback  and  three 
practical  controllers  for  a  more  complex  model  containing  six  bending 
modes.  The  results  are  not  markedly  different  using  a  more  complete 
model. 

Table  VI  shows  the  effectiveness  of  spoiler  controls  in  providing  gust 
load  relief.  These  are  for  state  controls.  Since  full-state  feedback  was 
used,  conservative  depreciation  was  used  from  these  answers  to  provide 
the  numbers  of  Table  I  representing  the  results  of  simplified  control. 

Table  vm  presents  results  for  maneuvering  load  control  (MLC).  The 
maneuver  is  a  step  column  input  that  attains  a  steady-state  value  of  1. 5 
incremental  g. 

The  free  aircraft  does  not  use  symmetric  ailerons  in  the  steady  state. 

At  +1. 5  Incremental  g,  the  wing  root  perturbation  stress  is  -17,400  psi. 

The  peak  perturbation  stress  achieved  during  the  transient  for  the  step 
column  input  is  -18,700  psi  (row  1  of  Table  Vin). 
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Controllers  3A  and  3B  are  based  on  the  Iteration  3  controller.  On  con* 
troller  3A,  af.taron  feedbacks  were  allowed,  with  significant  improvements 
in  stress  and  stress  rate  at  the  expense  of  bending- mode  damping. 

Controller  14  practical  controllers  14A,  14B,  and  14C  produced  lower 
stress  Sj  levels  than  the  others.  Stress  rate  levels  were  higher.  On  con¬ 
troller  14C,  measurement  6  (lagged  acceleration)  increased  bending-mode 
damping  and  lowered  the  short- period  frequency.  The  short-period  damping 
ratio  was  somewhat  higher,  however. 

Controllers  5D,  3B,  and  14C  are  the  most  desirable  controllers  from 
the  bending-mode  damping  viewpoint,  which  could  be  important  for  stability 
margins.  Controller  3B  is  the  least  complex  of  the  three)  however,  14  C  is 
only  different  by  a  lag  network.  Controller  5D  uses  an  extra  sensor. 

Table  V  summarizes  the  performance  of  two  state-feedback  and  three 
practical  controllers  for  a  more  complex  model  containing  six  bending 
modes.  The  results  are  not  markedly  different  using  a  more  complete 
model. 


Table  VI  shows  the  effectiveness  of  spoiler  controls  in  providing  gust 
load  relief.  These  are  for  state  controls.  Since  full-state  feedback  was 
used,  conservative  depreciation  was  used  from  these  answers  to  provide 
the  numbers  of  Table  I  representing  the  results  of  simplified  control. 

Table  vni  presents  results  for  maneuvering  load  control  (MLC).  The 
maneuver  is  a  step  column  input  that  attains  a  steady-state  value  of  1. 5 
incremental  g. 

The  free  aircraft  does  not  use  symmetric  ailerons  in  the  steady  state. 
At  +1. 5  incremental  g,  the  wing  root  perturbation  stress  is  -17,400  psi. 
The  peak  perturbation  stress  achieved  during  the  transient  for  the  step 
column  input  is  -18,700  psi  (row  1  of  Table  Vin), 


/ 
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Table  VIII.  Maneuver  Load  Control® 


*For  '•’1.9  incremental- 1  command.  Relief  la  linear  with  control. 
Hence,  reault*  may  be  ueed  to  determine  effecttveneee  with 
different  eurface  deflections. 


The  effects  of  connecting  the  Input  of  the  aileron  actuator  to  the  control 
column  are  shown  in  row  2  of  Table  VIII.  Steady- state  values  at  1. 5 
perturbation  g  of  Sj  and  5u  are  -9810  psi  and  25  deg.  Thus,  the  steady-state 
MLC  relief  is  1  -  ^ ■  0.4362.  The  ratio  of  peak  stresses  without  and 
with  MLC  for  the  frets  aircraft  is  1  -  ^ -  0.  3529. 

Table  VIII  (row  3)  shows  that  adding  the  spoiler  to  the  free  aircraft 
provides  further  reductions  in  both  steady  and  peak  wing  root  stresses. 

Table  VIII  (row  4)  shows  that  feedback  provides  attenuation  of  peak 
maneuvering  stresses;  but,  of  course,  it  can  do  no  better  than  the  free  air¬ 
craft  (with  MLC  ailerons)  in  the  steady  state  at  +1. 5  incremental  g. 
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SECTION  V 

LICS  FUNCTIONAL  BLOCK  DIAGRAMS 


The  block  diagram  corresponding  to  simplified  controller  14C  (Table  IV) 
is  shown  in  Figure  8.  The  MLC  feedforward  gain  was  computed  to  enforce 
the  appropriate  steady-state  aileron  surface  deflection  per  g.  At  a  single 
flight  condition  this  was  satisfactory.  However,  flight  throughout  the  enve¬ 
lope  would  have  required  excessively  complex  scheduling  because  of  the 
large  variation  in  control  column  deflection  per  g  with  center- of- gravity 
variations.  The  system  was  then  revised  to  provide  maneuvering  relief  with 
a  feedback  control  (Figure  9)  using  the  inputs  to  the  aileron  and  elevator 
servos.  The  only  feedforwards  remaining  are  the  existing  mechanical  links 
in  the  aircraft.  Gains  and  time  constants  were  not  determined  for  this  con¬ 
figuration. 

Similarly.  Figure  10  illustrates  a  possible  functional  block  diagram  of  a 
system  using  the  spoilers  in  addition  to  the  aileron  and  elevator.  The  various 
gains  shown  were  not  determined. 

System  I  (Figure  9)  requires  one  additional  pair  of  accelerometers  in 
the  wings  as  input  sensors,  and  it  makes  use  of  the  autopilot  normal  accel¬ 
erometers  and  the  pitch  augmentation  rate  sensors  already  on  board  the  air¬ 
craft.  System  II  (Figure  10)  requires  two  additional  pairs  of  dual  accelerom¬ 
eters  in  each  wing,  and  also  requires  the  addition  of  dual-hydraulic  servos  to 
control  the  spoilers.  A  norm-rough  weather  switch  is  also  required  in  the 
cockpit. 


Figure  8.  Simplified  Control  Block  Diagram 
(Controller  14C) 


Figure  10.  System  II  LICS  (With  Spoilers)  Block  Diagram 


SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


Results  of  the  LICS  study  showed  that  rms  wing  root  stress  and  stress 
rate  can  be  reduced  by  factors  of  SO  percent  and  31  percent,  respectively, 
by  using  symmetric  aileron  and  elevator  control.  However,  two  deficiencies 
become  apparent: 

#  Handling  qualities  are  degraded  somewhat  (this  is  noticeable 
when  observing  transients  from  step  commands). 

e  No  means  for  enforcing  steady* state  load  relief  are  included. 

At  the  time  this  report  was  being  written,  additional  work  had  been  done 
on  the  C-5A  under  contract  from  the  Air  Force. 

Results  from  this  study  (Ref.  6)  indicate  that  the  handling  qualities  can 
be  maintained  with  some  loss  in  rms  stress  performance.  The  second 
deficiency  can  be  corrected  by  at  least  two  techniques: 

e  Integral  control  can  be  used  to  enforce  steady-state  load 
relief  proportional  to  normal  acceleration. 

e  Direct  accelerometer- to* aileron  feedback  can  be  used  to  give 
steady-state  load  relief.  A  highpass  network  is  used  on  the 
aileron  input  to  wash  out  all  other  steady-state  inputs. 

;'l  ,  ■’]  i 

In  future  studies  it  is  recommended  that  wing  torsion  be  considered. 
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APPENDIX 

SYSTEM  MATRICES  FOR  FLIGHT  CONDITION  37 


For  the  usual  notation  the  model  has  the  form 
x  ■  Fx  +  GjU  +  Ggfl 
r  *  Hx  +  Du 
y  ■  Mx 

where 

xT  *  v  Tb*  ^1*  TV  *V  ®a*  °e  '  pl’  p2*  p3* 

p4,  P5,  P6.  Wg) 

T  #  •  a 

r  *  (Sj,  ®2 #  My  Mgy t  d/ nj ) 

and  y  is  defined  in  Table  VII. 

The  six  matrices  (F,  Gj,  G2#  H,  D,  and  M)  corresponding  to  flight  condition 

37  follow. 
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